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Abstract 

We consider a mass-asymmetric polarized Fermi system in the presence of Hartree-Fock (HF) 
potentials. We concentrate on the BCS regime with various interaction strengths and numerically 
obtain the allowed values of the chemical and HF potentials, as well as the mass ratio. The func- 
tional dependence of the heat conductivity of the N-SF interface on relevant physical parameters, 
namely the temperature, the mass ratio, and the interaction strength, is obtained. In particular, 
we show that the interface conductivity starts to drop with decreasing temperature at the tem- 
perature, T m , where the mean kinetic energy of the particles is just sufficient to overcome the SF 
gap. We obtain T m as a function of the mass ratio and the interaction strength. The variation of 
the heat conductivity, at fixed temperature, with the HF potentials and the imbalance chemical 
potential is also obtained. Finally, because the range of relevant temperatures increases for larger 
values of the mass ratio, we consider the 6 Li- 40 K mixture separately by taking the temperature 
dependence of the pair potential into account. 
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I. INTRODUCTION 



Recently, the study of the behavior of ultracold Fermi gases with two imbalanced hy- 
perfine states has opened up an interesting new area in many-body atomic physics. In this 
connection, extensive studies have been reported that propose various candidates for the 



pairing state, including the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state \i 1 
normal phase separation the Sarma state [4|, the p-wave pairing state 



deformed Fermi surface superffuid [7]. Of central importance is the phase separation of a 
superffuid (SF) paired core surrounded by a polarized normal (N) phase, where, in addition 



the BCS- 
|, |6j and the 



Such a phase-separation scenario had been proposed by Clogston 



to theoretical studies |8Ml9|. important experimental work has been carried out 



24| and Chandrasekhar 



25] long ago, who predicted the occurrence of a first-order transition from the N to the SF 



state. An interesting result in this connection is the appearance of a temperature difference 
between the two phases as a consequence of the blockage of energy transfer across the N-SF 
interface. This blockage is due to a SF gap, which causes low-energy normal particles to 
be reflected from the interface. By studying particle scattering off the interface, the heat 



conductivity has been calculated 



26 



In this paper we consider a polarized Fermi system consisting of two spin species with 
unequal masses in the presence of HF potentials. We concentrate on the BCS regime with 
various interaction strengths and numerically obtain the allowed values of the chemical and 
HF potentials, as well as the mass ratio m r . The functional dependence of the heat conduc- 
tivity of the N-SF interface on the relevant physical parameters, namely the temperature, 
the mass ratio, and the interaction strength, is studied in detail. Our focus is on energies 
slightly above the transmission threshold, because we are considering low temperatures. In 
our calculations, we therefore use the approximate low-temperature form of the Fermi-Dirac 
distribution and regard the pair potential as temperature-independent. In particular, we 
show that the interface conductivity starts to drop with decreasing temperature at the tem- 
perature, T m , where the mean kinetic energy of the particles is just sufficient to overcome the 
SF gap. The drop is, thus, a result of the blockage of the energy transfer due to the reflec- 
tion of particles from the interface and signifies a build-up of temperature difference across 
the interface. We obtain T m as a function of the mass ratio and the interaction strength. 
The variation of the heat conductivity, at fixed temperature, with the HF potentials and 
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the imbalance chemical potential is also obtained. Finally, we single out the particular case 
of the 6 Li- 40 K mixture (m r = 6.7), due to its importance in experimental and theoretical 



studies 129 



3l|. Because the range of relevant temperatures increases for larger values of the 
mass ratio, here we take the pair potential to be temperature-dependent and use the exact 
Fermi-Dirac distribution instead of its approximate low-temperature form. 

II. HARTREE-FOCK POTENTIAL AND MASS ASYMMETRY EFFECTS 

We consider a polarized Fermi gas consisting of two fermionic species (imbalanced hy- 
perfine states f, l) of masses m-\- , and chemical potentials fi^, ^4. at sufficiently low tem- 
perature. The "t — 4- interaction is assumed to be a contact interaction characterized by the 
coupling constant V = —4na/mR (h = 1) with niR = Im^m^j [m^ + mi). For the superffuid 
phase we define the species-imbalance chemical potential, h s = (/if — A*j.)/2, and the average 
chemical potential, ji s = (/i^- + //j.)/2 — U s , where U s is the superfluid HF potential. For 
the calculation of transmission coefficients, and hence the heat conductivity, we need to 
obtain the solution of Bogoliubov-de Gennes equations [32]. The effective hamiltonian of 
the system may be written as 



H = I d 3 xY^[^{ri)H(ri)il)(ri) + U(ri)4>\ri)ip(ri) 

i 

+A(r)V> t (r t)^ f (r |) + A*(r)^(r |)^(r f)] (1) 

where H{ri) = - ^ (i =t,|) and 

U(r |) = -V < ^(r l)ijj(r I) >, U(r |) = -V < ifl{r t)ip(r |) > 
A(r) = -V < ip(r i)ip(r f) >= V < ip(r f)^(r |) > 

are the HF and pair potentials, respectively. We use the approximation that U and A are 
independent of r. It is noted that in the superfluid phase (unlike the normal phase) all the 



HF potentials are equal 



33|. The traditional forms of ip(r t) and ip(r \) are 



t) = £(7*t«*(r t) - llvt(r |)), ^(r I) = 5>*«*(»" I) + 7 lt«J(r I)) (2) 
fc fc 



where 7,7' (u, v) are the fermionic quasiparticle operators (wavefunctions) . By using these 
expressions and the commutation relations between 7, 7^ and H, one can straightforwardly 
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obtain the Bogoliubov-de Gennes equations 

[H{r |) + U(t)]u(r f) + Av{r |) = Eu{r f) 
A*u(r |) - [H(r |) + U(\.)}v(r |) = Ev{r |). (3) 

One obtains the second set of equations by interchanging "\ and \.. In the a-channel, we take 
u(r f) for the particle- like and v(r \) for the hole-like wavefunctions. 

To proceed, let us take the N-SF interface to be in the x = plane and introduce the 
superscript s (n) for the solutions in the SF (N) phase. We also introduce <j>k( q ){ r ) = 
exp[z(fc|| • r ± kr q )x)] for the N phase, and 0^ (r) = exp[i(fcii ■ r ± fc^x)] for the SF phase, 
where q = p,h refers to particle, hole and k\\ (k^ q \ k( q )) denotes the component of the wave 
vector k parallel (perpendicular) to the interface. Notice that q appears as a subscript 
(superscript) for the N (SF) phase throughout our notation. 

From the Bogoliubov-de Gennes equations we obtain the relations, k? h \ = fcjL + 
2m t [C/(t) - m r U{i) + U s {m r - 1) - 2e) and k^ 2 = kf p) + 2m t (C/(t) - U a - f*), where 
m T = m^jm^ (mass ratio), 2e = (E + h 8 )(l + m r ) + /i s (l — m r ), and ^ = e =F V £2 ~ m r A 2 . 
Thus for the N phase we write 



ut\r t) = E WwW, 4 n V I) = £ WW r )" ( 4 ) 

<T=± <T=± 

As for the SF phase, 



4 S V t) = E UZ {9) tf 9 \r), vi s \r |) = E ^tfV) (5) 

where 1/ fc CT(p) = BU^ and (h) = B^uf^ with 5 = £+/A. The amplitudes etc. 
are to be determined by matching the wave functions and their derivatives at x — 0, of 



course 



34|. 



Denoting = k 2 p J2m^ for £ + -£/(t) + t/ s < £(p) < /u^ — Z7 (t)+-E, particle-like and hole- 
like excitations both occur in the SF side, but Andreev reflection [33] is forbidden. However, 
for m — C/(t) + E < ^ < 2e — Z7(t) + m r U(\) - U s (m r - 1), we have particle-like and 



hole-like excitations, as well as normal and Andreev reflections |27(. In other regions, the 
particle has insufficient energy to excite the SF side and, thus, the transmission coefficients 
vanish. We, therefore, restrict our attention to the above two regions, which we shall denote 
by I and II, respectively. Moreover, our focus is on energies slightly above the transmission 
threshold (e « y/rn^A), because we are considering low temperatures. 



Denoting the x-component of the current density by j x , the transmission coefficient is 
given by W = j^/jx, where the superscripts T and I refer to the transmitted and incident 
quasi-particle current densities, respectively. The general form of j (for a-channel) is 



3a{r) = -^W(r t)Vu(r t) - u(r t)Vu*(r t)] 

-v*(rl)Vv(rl)+v(rl)Vv*(rl)]. (6) 



2m± 

The heat conductivity (for a-channel) is given by 

K= n 2 (l + m r ) 2 df 1 1 d Z(p) d c(z-e )f(e,T)W(e^ {p) ) + (t^l,p^h) (7) 

where Eq = e\e=o and f(e,T) is the Fermi-Dirac distribution, which, in the low temperature 
limit T <C ^Jrf^.A (ks = 1), reduces to e _Tm//T (up to a proportionality constant), where 

-T m = 2^- -. 8 

2 m l + m r y J 

The right hand side is the minimal energy attained by the a spectra, which is positive (in 
order to have a gapped spectrum) and independent of the HF potentials. 

In the mass asymmetric case, analytical calculation of the heat conductivity in the BCS 
regime (unlike the deep BCS regime in which the Andreev approximation is valid) is a 
formidable task, especially when HF potentials are present. We, therefore, approach the 
problem numerically and examine the effect of the HF potentials and mass asymmetry in 
regions I and II. We begin by obtaining the allowed range of values for all the relevant 
parameters in the BCS regime. To this end, we use the following standard relations. The 
HF potential of the superfluid phase, obtained by using the fermionic anticommutation 
relations for 7 and 7T, is given by the number equation 

where (k = £k — Us = k 2 /2mn — fi s . Similarly, the gap equation is given by 

1 f d 3 k .11. . . 

1 = -V / - — - ; . 10 



The above integrals can be calculated using [36 ] 



00 »A 1 

dz - i = — !!— (1 + x 2 ) x ' 2 P x { , ) (11) 

[(z- l) 2 + x 2 ]V2 sinvrA^ ; y/T+lc 2 V ; 
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where Pa is the Legendre function. Equations (jHJ) and ({TO]) , thus, yield 

(12) 



2 _ ft/a to 



where = —[1 + (A//i s ) 2 ] x / 2 . Since n s = /cp/37r 2 , using ( !T2|) we find 

n ^~3 4 ^i 3 /2(0 / -1 q \ 

(k¥a) = ^ P^)-^nW (13) 

This relationship determines the allowed values of q by fixing 1/kpa in the BCS regime. 
Through ffT2l) we thus find U s , A, and the latter yields h s via the Clogston limit. We, 
therefore, have m and Hi as well. In the normal phase we similarly find [27J 

which yield the allowed values of U(f) and U(\.) in the BCS regime. For m r less than a 
cut-off value M (which depends on the interaction strength), we find two solutions for each 
U(i), satisfying U s < U(i) < (for m r > M, no real solution exists). Since the interactions 
are attractive the upper bound on the potentials is obvious. The lower bound implies that 
the density of the SF phase (the core region) exceeds that of the N phase. This reconciles 
with the fact that a N-SF interface exists, separating an unpolarized SF from a partially 
polarized N phase. Therefore, we take 1 < m r < M for the allowed range of values of m r in 
the BCS regime. 

III. RESULTS AND DISSCUSION 

The relevant physical parameters in a dilute Fermi gas are the temperature, the mass 
ratio, and the interaction strength. It is, therefore, important to find how physical quantities 



depend on these parameters [37|. For the heat conductivity, (0) is found to give 



^ = g (t) e 2T (15) 

where T m is defined by (jSJ), /% = T \ji^m^+ H\jn\) I 'vr 2 is the heat conductivity of the N phase, 
and Tp is the Fermi temperature. The functional forms of G(^,nv) and T m (^^,m r ) will 
be discussed shortly. We note that, as a consequence of ( TT5l) . k/kn starts to drop from 
its maximum value at T = T m with decreasing temperature. This signifies a build up of 
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FIG. 1: (Color online) T m versus m r . 



temperature difference across the interface, which be understood as follows. According to 
OH]), at T = T m , the mean kinetic energy of the particles is just sufficient to overcome the SF 
gap. We, therefore, expect a blockage of energy transfer at lower temperatures, resulting in 
the reflection of particles from the interface and hence a drop in the interface conductivity. 
T m is an increasing function of m r , which, for fixed m r , increases with the interaction strength 
too (Fig. d). 

Fig. [2] shows typical results for the temperature variation of k/kn using l/k^a = 
—0.84, —0.67. As seen, for fixed m r , the larger the absolute value of 1/kpa (i.e., the weaker 
the interaction), the larger is the heat conductivity. Also, in the presence of HF potentials, 
the maximum value of k is almost the same for all m r for sufficiently weak interactions. 

Furthermore, the k/k,-^ at fixed temperature decreases with m r (Fig. [3]), resulting in an 
increase in the temperature difference across the interface. This means that the characteristic 
relaxation time increases with m r . Note that for sufficiently high values of the mass ratio, 
k/kn is independent of the interaction strength, provided T/Tp < 0.03. 

Fig. H] shows typical curves of k/k-^ versus T at fixed h s , for m r = 1, 1.4. At fixed T, 
k/kn decreases with increasing h s , because of the growing interaction strength. The effect of 
h s (which can be controlled by the species population imbalance) on k is more pronounced 
in the presence of HF potentials, because the latter affect the threshold line (e = A^Jm T ) 
by changing e. However, as seen, the role of h s diminishes at sufficiently low temperatures. 

The functional forms of G and T m in (fl5l) have been determined by the method of least- 
squares fit. They are valid for the whole BCS regime and reproduce the above results very 
accurately: 
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FIG. 2: Interface conductivity versus temperature for 1/kpa = —0.84 (left) and -0.67 (right), with 
(top) and without (bottom) HF potentials. 
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FIG. 3: (Color online) Interface conductivity versus mass ratio at T/T-p = 0.05 and 0.03 (inset). 



3 1 
2 m fc 



, m r 



Tpm r 
1 + m r 



t (m r ) + - — ti(m r ) + (-r—) 2 t2(m r ) 
k-pa k-pa 



(16) 



where 



0.93 0.43 1.78 0.79 0.75 0.32 

go = -0.51 + , 9l = -1.00 + , g 2 = -0.44 + 



t = 0.16 + 



1.38 1.14 



\/m r 



-0.99 



1.35 0.03 
+ 



\/m r 



-0.34 + 



0.39 0.07 



\/m r 
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FIG. 4: (Color online) Interface conductivity versus temperature for m r = 1 (inset) and 1.4, with 
(left) and without (right) HF potentials. 




FIG. 5: (Color online) Dependence of interface conductivity on temperature and interaction 
strength for m r = 1.2 (top) and 1.4 (bottom). 

The resulting functional dependence of k/kn on the temperature and interaction strength 
is depicted in Fig. El 

The curves of k/kn versus HF potentials at constant temperature are shown in Fig. El As 
seen, the role of the potentials becomes more important as the temperature increases. Also, 
the heat conductivity decreases with U^/U^, which is understood because of the resulting 
increase in the scattering length (and, hence, cross section). 

As a subsidiary result, we may also point out that our numerical calculations show that 
the role of the incident particles (from the N side) with energies in region I (where An- 
dreev reflection does not occur) is much more important in k than that of the incident 
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FIG. 6: (Color online) Interface conductivity versus HF potentials for m r = 1.4. 



particles/holes with energies in region II. 



The mixture 6 Li- 40 K 



Here, we consider in more detail the particular case of the 6 Li- 40 K mixture (m r = 6.7), 
due to its importance in experimental and theoretical studies. Since T m increases with m r 
(Fig. [1]), for larger values of the mass ratio such as here, the range of relevant temperatures 
increases. Hence, it would be more appropriate to take the temperature dependence of A 
into account. We have 1381 



A 



1 oc (8- 



7\ 
A' 



-A/T 



(17) 



where A is the zero temperature limit considered in our previous calculations. Also, we take 
the distribution function f(e,T) in (J7J) to be the exact Fermi-Dirac distribution instead of 
its approximate low-temperature form. (However, for simplicity, we take HF potentials to 
be zero.) 

We find the following analytic expression for the transmission coefficient of region I: 
^ = 8(e - e )[2v^;A(r)(e - ^(T))]^ x 
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FIG. 7: (Color online) Interface conductivity versus temperature and interaction strength for 6 Li- 
40 K mixture. 
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FIG. 8: K m and r m versus interaction strength for 6 Li- 40 K mixture. 



\Jx(x-lW + ™?)(2 X - 3) + 2K - l)v/(x-l)(x-2)] 



(18) 



(m r - 1) + (Vx^Z ~ Vx)[( m r ~ 1) VX ~ (rnr + l)Vx=l}_ 

where x — £(p)/ A(T)y / m7- By using equations (fTTl) . (|T2l) . and the definition of <r, we can 
obtain the temperature dependence of fi s . The interface conductivity (j7j) is then obtained 
via numerical interpolation (Fig. [7]). Similarly, the functional forms of K m and r m (the 
maximum values of k/kn and T/Tp, respectively) have been determined by the method of 
least-squares fit. They are 



1 



kpa' 



0.006 + 0.02-?- + 0.02(-^) 2 
kpa kpa 

- 0.2 - 0.04^ 0.14(- 1 ^ 2 



(19) 



kpa k-pa kpa 

Fig. [8] shows the graphs of n m and r m versus the interaction strength. 

It is noteworthy that, since h s < A /rfvA(T) < fi s , where the second inequality sign is 
owing to the fact that k is real, the condition for Clogston limit (h s <C fi s ) is satisfied more 
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stringently as m r increases. 
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